Hawking radiation from ultrashort laser pulse filaments 
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Event horizons of astrophysical black holes and gravitational analogues have been predicted to 
excite the quantum vacuum and give rise to the emission of quanta, known as Hawking radiation. 
We experimentally create such a gravitational analogue using ultrashort laser pulse filaments and our 
measurements demonstrate a spontaneous emission of photons that confirms theoretical predictions. 



In 1974 Hawking predicted that the space-time curva- 
ture at the event horizon of a black hole is sufficient to 
excite photons out of the vacuum and induce a contin- 
uous flux, referred to as Hawking radiation [H, Q ■ In a 
simplified description of the process, vacuum fluctuation 
pairs close to the horizon are split so that the inner pho- 
ton falls in and the outside photon escapes away from 
the black hole. As the outgoing photon cannot return 
to the vacuum, it necessarily becomes a real entity, gain- 
ing energy at the expense of the black hole. It was soon 
realized that the essential ingredient of Hawking radia- 
tion was not the astrophysical black hole itself but rather 
the space-time curvature associated to the event hori- 
zon 043 ■ There are a wealth of physical systems that 
may exhibit event horizons ranging from flowing water or 
Bose-Einstein- Condensates to a moving refractive index 
perturbation (RIP) in a dielectric medium 0, 0, 
In a few words, referring to the case of optical pulses in a 
dielectric medium proposed by Philbin et al. |9j, a laser 
pulse with large intensity, /, propagating in a nonlinear 
Kerr medium will excite a RIP given by 5n = 712/ where 
tt-2 is the so-called nonlinear Kerr index [111 ]. Light ex- 



periences an increase in the local refractive index as it 
approaches the RIP and is thus slowed down. By choos- 
ing appropriate conditions (frequency of the light and 
velocity of the RIP) it is possible to bring the light waves 
to a standstill in the reference frame comoving with the 
RIP, thus forming a so-called white hole event horizon, 
i.e. a point beyond which light is unable to penetrate. 
A similar mechanism may be observed with water waves 
or with any kind of waves in a flowing medium, so that 
the formation of an analogue event horizon is a rather 
universal phenomenon that may be studied in accessible 
laboratory conditions. What remains to be established, 
is whether Hawking radiation is actually emitted in the 
presence of an event horizon of any kind, be it analogue 
or astrophysical. 

Recently, an alternative approach was proposed for gen- 



*Eloctronic address: E-mail: d.faccio@hw.ac.uk 



crating controllable RIPs, namely ultrashort laser pulse 
filamentation [Tfjj . 

Ultrashort laser pulse filaments are intense laser pulses in 
a transparent Kerr medium (i.e. with a third order opti- 
cal nonlinearity) characterized by a high-intensity spike 
that propagates apparently without diffraction over dis- 
tances much longer than the Rayleigh length associated 
to the spike dimensions and have been proposed for many 
applications ranging from white light g eneration to the 
control of atmospheric conditions |14| - tl8| . Filaments may 
either occur spontaneously when a powerful Gaussian 
shaped beam is loosely focused into the Kerr medium 
EH or, alternatively, they may be induced by pre- 
shaping the laser pulse into a Bessel beam [l9|, [2(| . 
In this Letter we use ultrashort laser pulse filaments to 
create a traveling RIP in a transparent dielectric medium 
(fused silica glass) and we report experimental evidence 
of photon emission that on the one hand bears the char- 
acteristics of Hawking radiation and on the other is dis- 
tinguishable and thus separate from other known pho- 
ton emission mechanisms. We therefore interpret the 
observed photon emission as an indication of Hawking 
radiation induced by the analogue event horizon. 
Neglecting any dependence on the transverse coordinates 
and dispersion it is possible to give a complete description 
of the event horizon associated to the RIP and calculate 
a blackbody temperature of the emitted photons in the 
laboratory reference frame [sHH)]. However, we point out 
that the dielectric medium in which the RIP is created 
will always be dominated by optical dispersion and this 
in turn implies that in real settings the spectrum will not 
be that of a perfect blackbody and that in any case only 
a limited spectral portion of the full spectrum will be 
observable. To show this last point we describe the RIP 
as a perturbation induced by the laser pulse on top of a 
uniform, dispersive background refractive index no, i.e. 
n(z,t,uj) = rio(w) + Snf(z — vt), where ui is the optical 
frequency, f(z — vt) is a function bounded by and 1, 
that describes the shape of the laser pulse. In the ref- 
erence frame co-moving at velocity v with the RIP, the 
event horizon in a 2D geometry is defined by c/v = n 
which admits solutions only for RIP velocities satisfying 
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FIG. 1: Experimental layout used for detecting analogue 
Hawking radiation. The input laser pulse is focused into a 
sample of fused silica (FS) using an axicon or lens (F). An 
imaging lens (I) collects the photons emitted at 90 deg and 
sends them to an imaging spectrometer coupled to a cooled 
CCD camera. 

the inequality [l2| 



no(u>) + Sn c no(uj) 

Therefore, Eq. ([T]) predicts an emission spectrum with 
well-defined boundaries and it is precisely this feature of 
the spectral emission that we consider to be peculiar to 
analogue Hawking radiation in the present settings. 

We note that two horizons may be associated to each 
RIP: the leading edge is the analogue of a black hole 
horizon whilst the trailing edge has time reversed fea- 
tures and is the analogue of a white hole horizon [j| ■ 
The experimental layout is shown in Fig. Q] The laser 
pulses are provided by a regeneratively amplified, 10 Hz 
repetition rate Nd:glass laser. The pulse duration is 1 ps 
and maximum energy is 6 mJ. The Bessel pulse filament 
with a cone angle = 7 deg is generated by a 20 degree 
fused silica axicon (conical lens), F, placed directly in 
front of the 2 cm fused silica Kerr sample where the RIP 
is generated. The input energy is varied in the 100-1200 
//J range. 

Radiation from the filament is then collected at 90 de- 
grees with respect to the laser pulse propagation axis us- 
ing a lens that images the filament on to the input slit of 
an imaging spectrometer. The spectrum is then recorded 
with a 16 bit, cooled CCD camera. This arrangement, in 
particular detection at 90 degrees, was chosen in order to 
strongly suppress or eliminate any spurious effects. More 
specifically: 

(i) Cerenkov-like radiation, i.e. in this context, radiation 
from a superluminal perturbation: spontaneous emission 
from a strictly superluminal perturbation has been con- 
sidered in detail in Ref. [Hj]. The most relevant differ- 
ence with respect to the present measurements is that 
Cerenkov-like emission occurs with no upper bound in 
the spectral emission window, in contrast to the limited 
Hawing spectrum given by Eq. (JXJ) . 

(ii) Four wave mixing (FWM) and self phase modula- 
tion (SPM) will not occur at 90 degrees due to phase- 
matching constraints which imply that any newly gener- 
ated frequencies will be generated at small angles with re- 
spect to the propagation axis [2l| . In any case we directly 



verified that for the relatively large = 7 deg Bessel cone 
angle used in our experiments, no FWM, SPM or spec- 
tral broadening was observed at any angle, even in the 
forward direction. 

(iii) Rayleigh scattering will occur at 90 degrees only for 
vertically polarized light and the scattering process will 
maintain the polarization state. In our experiments we 
used horizontally polarized light and in any case, in virtue 
of point (i), there is no input or generated light at the 
frequencies relevant for the present experiment. 

(iv) Fluorescence is certainly the main problem in these 
experiments. It bears many features in common with 
Hawking radiation yet, it may still be clearly distin- 
guished from the latter. Figure HJa) shows the overall 
spectrum measured at 90 degrees and integrated over 30 
laser pulse shots. In the figure, R, Fl and F2 indicate 
the laser pulse induced spontaneous Raman, nonbridg- 
ing oxygen hole centre (NBOHC) and oxygen deficient 
centre (ODC) fluorescences, respectively. These fluores- 
cence peaks are well documented features of fused sil- 
ica [12, HH . This therefore allows to fit the peaks with 
Gaussian functions (in the frequency domain) and subse- 
quently subtract out the fluorescence signals, thus lead- 
ing to greatly improved contrast and cleaner spectra. 
Nevertheless, looking at the fluorescence spectrum it is 
clear that, given the possibility to tune the Hawking emis- 
sion window by tuning the group velocity of the laser 
pulse, the most advantageous situation is that in which 
the window is located between 800 nm and 900 nm where 
there is no known fluorescence emission and moreover the 
CCD response is maximum. In fused silica, n<x ~ 3-10~ 16 
cm 2 /W and typical intensities of / ~ 10 12 — 10 13 W/cm 2 
may be obtained, thus leading to a Sn ~ 10~ 2 — 10 ~ 4 . 
Therefore, if we take for example a Gaussian pulse cen- 
tred at 1055 nm, with Sn = 0.001 and group velocity, 
vq = dco/dk, determined solely by material dispersion, 
then according to Eq. ([I} wc would expect emission in 
the region between 500 nm and 510 nm. This would 
clearly fall in between the Fl and F2 fluorescence peaks 
and leads to noisy results even after subtraction of the flu- 
orescence signal. By using Bessel pulses we were able to 
solve this problem. Figure [^b) shows the Hawking spec- 
tral window predicted from Eq. ([1]), where the dashed 
line shows the value c/vb with = 7 deg and the Bessel 
pulse velocity is vb = vgJ cos 0. This line intersects the 
refractive index curves n and n + Sn with Sn = 0.001 
delimiting a window (shaded area) which now lies in 
the desired region, i.e. between 800 nm and 875 nm. 
Figure (Ha) shows the resulting spectra, integrated over 
3600 laser shots. The black line shows a reference spec- 
trum obtained from a Gaussian pulse, that clearly shows 
the absence of any signal notwithstanding the same peak 
intensity of the Bessel pulses and thus underlining the 
absence of any fluorescence or other possible emission 
signals when the RIP v does not satisfy Eq. (Q}. The 
other four curves show the emitted spectra for four dif- 
ferent input energies of the Bessel filament, as indicated 
in the figure. A clear photon emission is registered in the 



3 



600 

c 

5 300 



- 1 1 1 1 


i i i i 


1.457; 




1 (a) 




Fl 










1.453^ 


- F2 J 










1.449- 





200 600 1000 
wavelength (nm) 




700 900 1100 
wavelength (nm) 



FIG. 2: (a) Measured CCD photo-electron (p-e) counts gener- 
ated by the fused silica fluorescence spectrum, (b) Prediction 
of the Hawking emission spectral range for the case of fused 
silica: no (blue solid line) and no + Sn (green dashed line) with 
Sn = 0.001. The shaded area indicates the spectral emission 
region predicted for a Bessel filament. 



wavelength window predicted by Eq. ([I]). Moreover we 
verified that the emitted radiation was unpolarizcd (data 
not shown) thus further supporting the interpretation of 
a spontaneous emission. In addition, the measurements 
clearly show that the bandwidth of the emission is in- 
creasing with the input energy. The Bessel pulse inten- 
sity evolution along the propagation direction z may be 
estimated analytically from the input energy using the 
2irIokzt&n 2 (8) exp 



-(z 2 tan 2 0)/w 2 ] 



24] 



formula I(z) 

where k — uj/c and I and \/2wo are the input Gaussian 
peak intensity and radius at 1/e 2 , respectively. By fitting 
the measured spectra with Gaussian functions [dashed 
curves in Fig. [3f a)] we may therefore estimate the band- 
width as a function of input energy and Bessel inten- 
sity. Moreover, using the fused silica dispersion relation 
shown in Fig. UtR) it is possible to map the bandwidth 
into values of Sn so that we finally obtain Fig. [3jb) that 
shows the bandwidth and the Sn as a function of input 
energy and Bessel pulse peak intensity (at z = 1 cm 
where measurements were performed). There is a clear 
linear dependence which is in qualitative agreement with 
the fact that, according to our interpretation, the emis- 
sion bandwidth is predicted to depend on Sn which in 
turn is a linear function Sn = nil of the pulse intensity. 
The slope of the linear fit gives n 2 = 2.8 ± 0.5 x 10~ 16 
cm 2 /W which is in good agreement with the tabulated 
value of n 2 ~ 3 x 10" 16 cm 2 /W [HEH. Therefore this 
shows that there is also an agreement at the quantitative 
level between the measurements and the model based on 
Hawking-like radiation emission. 

In Fig. Uwe present additional data regarding Hawking 
emission from a spontaneous filament obtained by replac- 
ing the axicon with a 20 cm focal length lens and thus 
loosely focusing a 50 uJ Gaussian pulse into the fused 
silica sample. The spontaneous nonlinear Kerr dynamics 
lead to the formation of a filament [l4[ that is very similar 
to the Bessel filament in the sense that it is characterized 
by a very high-intensity peak that propagates over long 
distances, thus creating a strong RIP in the medium that 
may be expected to excite Hawking radiation in a simi- 
lar fashion to that illustrated above. We underline some 
important differences with respect to the Bessel filament: 
the spontaneous filament is characterized by an intense 
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FIG. 3: (a) Spectra generated by a Bessel filament. Five 
different curves are shown corresponding to a reference spec- 
trum obtained with a Gaussian pulse (black line) and indi- 
cated Bessel energies (in fj,J) . Dashed lines are guides for the 
eye. The solid line connects the spectral peaks, highlighting 
the ~ 40 nm shift with increasing energy, in close agreement 
with the predicted 45 nm shift, (b) Bandwidths at FWHM 
and RIP Sn versus input energy and intensity. Solid line: 
linear fit Sn = nil with n 2 = 2.8 x 10~ 16 cm 2 /W. 



pulse that travels slower than the input Gaussian pulse, 
i.e. v <vq (sec 26] for a detailed description). Moreover 
the pulse accelerates during propagation so that a certain 
range of velocity values is covered during propagation, 
with low velocity at the beginning of the filament and 
higher velocity at the end. This has the effect of broaden- 
ing the emission window, which is predicted from Eq. ([1]) 
to be between 270 nm and 450 nm. The light gray-curve 
in Fig. Ufa) shows the spectrum measured when collect- 
ing light from the whole filament. As can be seen, the 
measured spectral bandwidth agrees well with the pre- 
diction. We also note that no emission was observed in 
this region with the Bessel filament, indicating that this 
is not a fluorescence signal. Moreover, by shifting and 
closing the input slit of the imaging spectrometer, we 
could collect spectra of the ending and beginning sec- 
tions of the filament, that are shown as dark-filled curves 
in Figs.fJJa) and (b), respectively. Remarkably, different 
sections of the filament emit only at selected portions of 
the overall spectrum. We underline that, accounting for 
the different pulse velocities in the different filament sec- 
tions, this striking behavior is in quantitative agreement 
with the predictions of Eq. (p}. 

Finally, we note that Eq. ([1]) and subsequent analy- 
sis, refers to the existence of phase velocity horizons. 
However, other analogue systems rely on the existence 
of group velocity horizons, i.e. points at which v = vq 
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FIG. 4: Spectra generated by the spontaneous filament, (a)- 
(b) shaded areas: spectra measured for two different positions 
of the imaging spectrometer input slit. The light gray curves 
in both figures show the spectrum measured with the input 
split fully open. The filament, imaged from the side at 90 
deg., is shown in (c) and (d). The vertical white lines show 
the position of the input slit. 



(see e.g. [23, [28|) so that there is still an open question 
regarding the relative role of the two different types of 
horizon. In the specific case of the Bessel pulse used in 
our experiments, a group horizon does not even exist, i.e. 
there are no frequencies that satisfy vg{u) = vb- Out 
measurements therefore indicate that, at least in the RIP 
setting, the phase velocity horizon alone may lead to pho- 
ton emission in the predicted spectral window. 
In summary, the same physics that underlie black hole 
evaporation in the form of Hawking radiation may be 
found and studied in other, more accessible systems. Our 
measurements highlight spontaneous emission of Hawk- 
ing radiation from an analogue event horizon generated 
by an "evaporating" refractive index perturbation and 
suggest a path towards the experimental study of phe- 
nomena traditionally relegated to the areas of quantum 
gravity and astrophysics. 
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